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The  mechanism  of  H20  dissociation  as  well  as  the  adsorption  and  oxidation  reaction  of  H2  on  yttria- 
stabilized  zirconia  (YSZ),  commonly  used  as  part  of  solid  oxide  fuel  cell  (SOFC)  anodes,  was  investigated 
employing  temperature-programmed  desorption  (TPD)  spectroscopy  and  density  functional  theory 
(DFT).  In  agreement  with  theory  the  experimental  results  show  that  interaction  of  gaseous  H20  with 
YSZ  results  in  dissociative  adsorption  leading  to  strongly  bound  OH  surface  species.  In  the  interaction 
of  gaseous  H2  with  an  oxygen-enriched  YSZ  surface  (YSZ  +  O)  similar  OH  surface  species  are  formed  as 
reaction  intermediates  in  the  H2  oxidation.  Our  experiments  showed  that  in  both  the  H20/YSZ  and  the 
H2 /YSZ  +  O  heterogeneous  reaction  systems  noticeable  amounts  of  H20  are  “dissolved”  in  the  bulk  as 
interstitial  hydrogen  and  hydroxyl  species.  The  experimental  H20  desorption  data  is  used  to  access  the 
accuracy  of  the  H2/H20/YSZ  adsorption/desorption  and  surface  reaction  kinetics  data,  employed  in  previ¬ 
ous  modeling  studies  of  the  electrochemical  H2  oxidation  on  Ni-pattern/YSZ  model  anodes  by  Vogler  et  al. 
[J.  Electrochem.  Soc.,  156  (2009)  B663]  and  Goodwin  et  al.  [J.  Electrochem.  Soc.,  156  (2009)  B1004].  Finally 
a  refined  experimentally  validated  H2/H20/YSZ  adsorption/desorption  and  surface  reaction  kinetics  data 
set  is  presented. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Among  the  different  types  of  fuel  cells,  all  of  which  convert 
chemical  energy  of  a  fuel  gas  directly  into  electrical  energy,  solid 
oxide  fuel  cells  (SOFCs)  possess  the  highest  electrical  efficiency. 
Furthermore  for  SOFCs,  as  compared  to  conventional  fuel  combus¬ 
tion  concepts  such  as  internal  combustion  engines  and  turbines, 
the  emission  of  e.g.  nitrogen  oxides  (NO*)  is  only  marginal.  Hence, 
SOFCs  are  widely  regarded  as  environmentally  friendly  energy 
sources  (see  e.g.  [1,2]  and  references  therein).  In  addition,  if  oper¬ 
ated  in  the  reverse  mode  as  solid  oxide  electrolysis  cells  (SOECs) 
they  can  be  used  to  store  electrical  power  at  times  of  electricity 
over-plus  [3].  However,  although  the  traditional  nickel  (Ni)/yttria- 
stabilized  zirconia  (YSZ)  anode  material  combination  is  widely  used 
in  SOFCs,  the  actual  microscopic  details  of  the  electrochemical 
reaction  mechanism  occurring  at  the  three-phase  boundary  (TPB) 
of  gas-phase,  electrode,  and  electrolyte,  are  not  yet  fully  under¬ 
stood  even  for  systems  as  simple  as  the  electrochemical  oxidation 
of  hydrogen  (H2),  which  can  be  described  by  the  global  reaction 
H2(g) +  02-ysz  H20(g)  +  2e-Ni,  and  the  electrochemical  reduction 
of  water,  H20(g)  +  2e_Ni  ->  H2(g)  +  02-Ysz>  respectively.  The  latter 
one  being  important  if  the  cell  is  operated  as  SOEC. 
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The  elucidation  of  the  surface  reaction  mechanisms  underly¬ 
ing  these  two  global  charge  transfer  reactions  is  considered  as 
one  of  the  most  difficult  aspects  of  understanding  the  fundamen¬ 
tal  electrochemical  behavior  of  SOFCs  [2].  As  outlined  in  a  recent 
review  article  by  Sun  and  Stimming  [4]  a  detailed  understanding 
of  the  anode  charge  transfer  reaction  mechanisms  and  kinetics 
plays  a  vital  role  in  the  systematic  development  and  optimiza¬ 
tion  of  advanced  SOFC  anode  materials.  In  recent  years  numerical 
simulations  based  on  computational  reaction-diffusion  models, 
which  explicitly  consider  the  coupled  behavior  of  reagent/product 
adsorption/desorption  and  diffusion,  elementary  thermal  surface 
reactions  and  elementary  charge  transfer  reaction  steps,  have 
emerged  as  valuable  tools  to  shed  light  on  the  electrochemical  reac¬ 
tion  processes  even  for  complex  geometries  of  ceramic-metallic 
composite  anodes  generally  employed  for  SOFCs  [5].  If  the  aim  is 
on  the  determination  of  quantitative  elementary  reaction  kinetics 
information  about  the  TPB  electrochemistry,  however,  numeri¬ 
cal  simulations  of  electrochemical  characterization  measurements 
(such  as  polarization  and  impedance  data)  of  model  anode  sys¬ 
tems  with  simplified  and  well  characterized  TPB  geometries  - 
such  as  patterned  and  point  electrodes  -  are  required.  At  this 
point  it  should  be  noted,  however,  that  the  reliability  of  the 
charge  transfer  reaction  kinetics  data  derived  with  this  approach 
strongly  depends  on  the  employed  electrode  kinetics  model  and 
the  accuracy  of  the  respective  elementary  surface  reaction  kinetics 
data  set. 
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Various  H2/H20/Ni-YSZ  electrode  kinetics  models  have  been 
reviewed  in  [4].  In  addition,  a  critical  summary  of  the  H2/H20/Ni- 
YSZ  electrode  kinetics  models  as  used  in  different  kinds  of 
numerical  simulation  approaches  of  experimental  SOFC  anode  data 
can  be  found  in  [6].  As  outlined  in  [4],  the  available  electrode  kinet¬ 
ics  models  mainly  disagree  with  regard  to  the  actual  location  where 
-  either  mainly  on  the  Ni  [7]  or  both  on  the  Ni  and  on  the  YSZ  surface 
[8,9]  -  the  heterogeneous  chemical  and  electrochemical  elemen¬ 
tary  reaction  steps  take  place.  As  noted  in  [6],  the  contribution  of 
the  reaction  of  H2  with  O2- Ysz  directly  on  the  YSZ  surface  -  followed 
by  electron  migration  to  Ni  through  the  YSZ  -  was  suggested  in  Refs. 
[10,1 1  ].  The  formation  of  interstitial  hydrogen  and  hydroxyl  at  the 
Ni  electrode/YSZ  electrolyte  two  phase  boundary  was  proposed  in 
Refs.  [12,13].  The  partial  contribution  of  water  production  and  dis¬ 
sociation  through  the  formation  of  OH  surface  species  on  YSZ  was 
put  forward  in  Refs.  [8,9].  Also,  as  have  been  shown  by  Sakai  et  al. 
the  effect  of  water  dissolution  in  ceria-zirconia-yttria  polycrystals 
may  influence  significantly  the  surface  exchange  rate  constant  [14]. 

Based  on  experimental  electrochemical  data  obtained  for  well- 
defined  Ni-pattern/YSZ  model  anodes,  which  revealed  that  the 
presence  of  H20  in  the  fuel  gas  has  a  “catalytic”  effect  on  the  anode 
kinetics,  a  mechanism  (similar  to  the  one  previously  proposed  in 
Ref.  [15])  was  suggested  where  the  presence  of  OH  on  the  YSZ 
surface  would  allow  for  an  efficient  proton  migration  on  the  YSZ 
(towards  the  Ni  electrode)  resulting  in  an  effective  enlargement 
of  the  electrochemical  active  region  for  the  anode  kinetics  [16].  In 
a  “Space-state  modeling”  analysis  of  the  latter  experimental  data, 
however,  due  to  the  lack  of  quantitative  adsorption  and  reaction 
kinetics  data  for  H  and  OH  on  YSZ,  the  authors  employed  a  mecha¬ 
nism  where  the  heterogeneous  reactions  involving  H,  OH,  and  H20 
take  place  solely  on  the  Ni  surface  [17].  As  a  consequence,  how¬ 
ever,  the  proposed  YSZ  mediated  catalytic  effect  of  H20  could  not 
be  investigated  in  the  framework  of  this  modeling  work. 

In  a  subsequent  more  detailed  analysis  of  the  Ni-pattern/YSZ 
model  anode  data  of  Ref.  [16],  which  was  based  on  an  elemen¬ 
tary  kinetics  surface  reaction-diffusion  model,  the  adsorption  and 
dissociation  reactions  of  H20  on  the  YSZ  surface  were  included 
[18].  A  sensitivity  analysis  performed  in  the  course  of  the  latter 
study  clearly  emphasized  the  great  influence  of  the  thermodynamic 
and  reaction  kinetics  data  of  YSZ  surface  species  on  the  electro¬ 
chemical  modeling  results  (see  e.g.  Fig.  5  of  Ref.  [18]).  A  similar 
surface  reaction-diffusion  model  based  computational  approach 
was  applied  in  Ref.  [19]  to  analyze  Ni-pattern/YSZ  model  anode 
data  reported  in  Refs.  [20,21].  The  surface  reaction  mechanism  of 
Ref.  [19]  included  the  adsorption  and  dissociation  reactions  of  H20 
as  well  as  H20  formation  due  to  H2  oxidation  on  the  YSZ  electrolyte 
surface.  In  both  computational  studies  [18,19],  good  agreement 
with  the  respective  experimental  target  data  [16,20,21],  could  be 
achieved,  although  the  H20/YSZ  adsorption/dissociation  and  sur¬ 
face  reaction  kinetics  data  were  considerably  different.  While  the 
data  employed  in  [18]  was  based  on  estimates  partially  guided  by 
results  of  quantum  chemical  calculations  of  H20  dissociation  on 
Ti02  [22],  the  kinetics  data  of  [19]  was  derived  by  the  authors  from 
a  sophisticated  computational  modeling  of  a  thermogravimetrical 
data  set  reported  in  Fig.  1 7  of  Ref.  [  1 5  ] .  In  the  light  of  the  importance 
the  H2/H20/YSZ  surface  kinetics  data  play  in  the  reaction  mecha¬ 
nisms  employed  in  the  computational  studies  of  Ni-pattern/YSZ 
model  anode  data  [18,19],  further  studies  are  desirable,  the  results 
of  which  allow  for  a  validation  and,  if  necessary,  the  refinement 
of  the  employed  H2/H20/YSZ  kinetics  data  sets.  Experimentally 
validated  kinetics  data  would  also  allow  to  assess  in  more  detail 
both  the  accelerating  effect  of  H20  on  SOFC  performance  as  well 
as  the  validity  of  the  assumption  underlying  the  computational 
approaches,  namely  that  the  charge  transfer  reaction(s)  are  the 
main  rate  limiting  steps  (for  further  details  see  [6]  and  references 
therein). 


Fig.  1.  Schematic  illustration  of  the  background  subtraction  procedure  applied  to 
account  for  H20  MS  signal  contributions  due  to  H20  desorbing  from  the  sample 
holder.  Curve  a:  Raw  H20  TPD  spectrum.  Curve  b:  Background  contribution  from 
the  sample  holder.  Curve  c:  Background  subtracted  TPD  spectrum  representing  H20 
desorbing  from  the  YSZ  sample.  The  TPD  spectra  are  shifted  vertically  against  each 
other  for  ease  of  clarity. 

In  the  following  results  of  temperature-programmed  desorp¬ 
tion  (TPD)  of  H20  dissociation/desorption  and  H2  oxidation  over 
YSZ  will  be  presented,  which  allow  for  a  direction  comparison  with 
simulated  TPD  spectra  based  on  the  H2/H20/YSZ  surface  kinet¬ 
ics  data  sets  of  Refs.  [18,19],  respectively.  In  addition,  following 
the  approach  of  Ref.  [23],  quantum  chemical  calculations  based  on 
density  functional  theory  (DFT)  were  performed,  to  investigate  the 
molecular  energetics  of  H20  adsorption/dissociation  and  H2  oxida¬ 
tion  on  YSZ. 

2.  Experimental  methodology 

The  TPD  studies  of  H20  desorption  and  H2  oxidation  on  YSZ 
were  performed  in  a  vacuum  chamber  (base  pressure  10-9  mbar), 
which  has  been  used  previously  to  investigate  CO  desorption  and 
oxidation  on  Pt,  Rh,  Ni  and  YSZ  surfaces  (for  further  details  see 
[24]  and  references  therein).  The  chamber  was  pumped  by  a  turbo 
pump  backed  by  a  Roots  pumping  station.  It  was  equipped  with  a 
quadrupole  mass  spectrometer  (QMS)  for  TPD  measurements,  and  a 
retarding  field  analyzer  for  low-energy-electron  diffraction  (LEED) 
studies.  In  the  present  work  polycrystalline  YSZ  sample  plates  (8.5% 
mol  Y203)  with  a  thickness  of  0.2  mm  were  employed.  The  origi¬ 
nal  YSZ  sample  plate  was  obtained  from  Itochu  (Tokyo,  Japan).  The 
temperature  of  the  YSZ  samples  could  be  varied  from  300  K  up  to 
1200  K.  The  substrate  temperature  was  measured  by  a  thermocou¬ 
ple  attached  to  the  sample  surface. 

Different  kinds  of  H20  TPD  measurements  were  performed 
in  the  course  of  the  present  work.  In  a  first  series  of  H20 
TPD  experiments  YSZ  samples  exposed  to  ambient  air  humidity 
(p(H20)  ^  24  mbar)  were  introduced  into  the  vacuum  chamber  and 
flashed  to  a  temperature  of  700  K  in  order  to  desorb  impurities 
and  physisorbed  H20  from  the  YSZ  sample  surface.  After  cool¬ 
ing  down  to  300  K  the  actual  desorption  spectrum  was  obtained 
by  recording  the  H20  mass  spectrometer  (MS)  signal  while  the 
sample  was  heated  from  300  K  to  ca.  1200K  with  a  linear  heat¬ 
ing  rate  of  1  Ks-1.  A  typical  raw  experimental  TPD  spectrum  is 
shown  in  the  upper  part  of  Fig.  1  as  curve  a.  In  order  to  obtain  a 
TPD  spectrum  representing  the  actual  desorption  of  H20  from  the 
YSZ  sample  the  H20  signal  background  originating  from  sample 
holder  parts  (curve  b  in  Fig.  1 ),  which  are  heated  up  with  some  time 
delay,  had  to  be  subtracted.  The  background  corrected  TPD  spec- 
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trum  originating  solely  from  H2O  molecules  desorbing  from  the  YSZ 
sample  is  depicted  as  curve  c  in  Fig.  1 .  All  TPD  spectra  shown  in  the 
following  were  obtained  by  applying  this  background  correction 
procedure. 

3.  Quantum  chemical  calculations 

Quantum  chemical  calculations  based  on  density  functional  the¬ 
ory  (DFT)  were  performed  to  investigate  H2  oxidation  and  H20 
dissociation  energetics  on  YSZ  surfaces.  As  in  our  previous  DFT 
studies  [25]  the  present  calculations  were  performed  using  the 
CASTEP  (Cambridge  Sequential  Total  Energy  Package)  computer 
code  [26]  in  the  framework  of  the  generalized  gradient  approxi¬ 
mation  (GGA),  as  proposed  by  Perdew  et  al.  [27],  in  combination 
with  Vanderbilt  ultrasoft  pseudopotentials  [28].  For  systems  with 
an  even  number  of  electrons  restricted  nonspin-polarized  cal¬ 
culations  were  performed.  For  systems  with  an  odd  number  of 
electrons  unrestricted  spin-polarized  calculations  were  performed. 
The  plane  wave  basis  set  was  truncated  at  a  kinetic  energy  of 
420  eV.  Computations  were  performed  over  a  range  of  k-points 
within  the  Brillouin  zone  as  generated  by  the  full  Monkhorst-Pack 
scheme  [29]  with  a  2x2x1  mesh.  A  further  increase  of  the  cut 
off  energy  and  the  number  of  k-points  resulted  in  negligibly  small 
changes  in  adsorption  energies  (typically  less  than  5kJmol_1), 
indicating  that  the  energy  values  are  well  converged.  A  slab  was 
repeated  periodically,  leaving  at  least  a  10  A  wide  vacuum  region 
in  the  direction  perpendicular  to  the  surface.  The  atoms  of  the  bot¬ 
tom  multilayer  were  fixed  in  all  calculations,  while  the  positions 
of  atoms  of  the  two  outermost  multilayers  and  the  positions  of 
atoms  of  the  adsorbed  molecules  were  optimized  during  the  pro¬ 
cess  of  structural  relaxation.  A  complete  linear  synchronous  transit 
(LST)/quadratic  synchronous  transit  (QST)  scheme  [30]  was  used  in 
the  transition  state  searches  for  the  surface  reaction  steps.  A  force 
convergence  criterion  of  0.03  eV  A-1  was  applied  in  all  calculations. 

The  slabs  used  in  the  present  study  to  investigate  H2  oxidation 
on  an  oxygen-enriched  YSZ  surface,  denoted  as  YSZ  +  O  in  the  fol¬ 
lowing,  and  H20  dissociation  on  a  YSZ  surface  are  reproduced  in 
Fig.  2(a)  and  (b),  respectively.  The  procedure  applied  to  obtain  a 
representative  slab  for  the  (111)  surface  of  YSZ  as  shown  in  Fig.  2(b) 
is  similar  to  the  one  described  in  Ref.  [23],  which  starts  from  a  pure 
(1 1 1)  terminated  zirconia  slab,  which  consists  of  12  Zr02  formula 
units  with  36  atoms  (for  further  details  see  Ref.  [23]  and  Fig.  2 
therein).  From  the  pure  zirconia  slab  an  oxygen  atom  was  removed 
to  account  for  vacancy  formation  (denoted  as  “vac.”  in  Fig.  2(b)) 
and  two  zirconium  atoms  were  substituted  by  two  yttrium  atoms 
which  are  depicted  as  black  atoms  (Y)  corresponding  to  a  9%  mol 
concentration  of  yttria  in  YSZ.  The  energetically  most  stable  YSZ 
slab  corresponds  to  the  structure  with  two  atoms  of  yttrium  close 
to  each  other  and  preferentially  occupying  the  next-nearest  neigh¬ 
bor  sites  to  the  oxygen  vacancy  [31,32].  Therefore,  for  calculations, 
we  adopted  the  mentioned  YSZ  slab  with  the  subsurface  vacancy 
and  two  Y  atoms,  substituting  Zr  in  the  outermost  and  the  sec¬ 
ond  surface  layers.  Finally,  in  order  to  generate  an  oxygen-enriched 
surface,  an  additional  oxygen  atom  (marked  by  *  in  Fig.  2(a))  was 
introduced  to  occupy  the  subsurface  oxygen  vacancy  of  YSZ.  The 
four  different  surface  oxygen  atoms  of  the  YSZ  +  O  surface  unit  cell 
are  indicated  by  the  numbering  (i  =  1-4). 

In  Fig.  3  two  low  energy  reaction  pathways  for  the  heteroge¬ 
neous  H2  oxidation  reaction  on  the  YSZ  +  O  surface  as  obtained  in 
the  DFT  calculations  are  schematically  depicted.  Both,  the  pathway 
shown  in  the  upper  part  (I)  and  the  pathway  shown  in  the  lower  part 
(II)  of  Fig.  3  lead  to  direct  formation  of  gaseous  FI20(g)  and  YSZ  by 
barrier  less  dissociative  adsorption  of  FI2  on  surface  oxygen  atoms 
via  the  formation  of  OH  surface  intermediates,  which  can  further 
react  to  yield  surface  H20  that  finally  desorbs  into  the  gas  phase 
leaving  back  an  YSZ  surface.  The  energetics  of  the  two  reaction 


Fig.  2.  Slabs  employed  in  the  present  DFT  calculations:  (a)  Oxygen  rich  YSZ  +  O  with 
the  additional  oxygen  atom  (*)  occupying  the  subsurface  oxygen  vacancy  of  YSZ,  the 
structure  of  which  is  reproduced  in  (b)  with  the  positing  of  the  subsurface  oxygen 
vacancy  marked  as  “vac.”  For  further  details  see  text. 

pathways  I  and  II  illustrated  in  Fig.  3  is  depicted  in  Fig.  4(a)  and  (b), 
respectively.  The  latter  results  also  reveal  that  the  reverse  reaction 
namely  the  dissociative  adsorption  of  H20  on  YSZ  leading  to  sur¬ 
face  OH  proceeds  is  an  activated  process  with  overall  DFT  activation 
energies  of  £a(I)  =  1 21 .3  kj  mol”1  and  £a(II)  =  91 .4  kj  moH1 .  The  lat¬ 
ter  value  is  in  reasonable  agreement  with  the  value  of  86.8  kj  mol-1 
in  the  sticking  probably  expression  used  by  Goodwin  et  al.  [19] 
to  describe  the  overall  process  of  dissociative  H20  chemisorption. 
In  the  mechanism  of  Ref.  [18],  on  the  other  hand,  the  dissociative 
chemisorption  of  H20  is  assumed  to  occur  with  considerably  lower 
activation  energy  of  9.6  kj  mol-1 . 

At  this  point  we  would  like  to  note  that  we  also  performed  com¬ 
prehensive  DFT  calculations  to  investigate  the  adsorption  of  H2  on 
the  YSZ  surface.  For  this  reaction  system,  however,  we  were  unable 
to  locate  any  low  energy  dissociative  adsorption  pathway  result¬ 
ing  in  surface  OH  formation  without  applying  artificial  constrains. 
Our  unconstrained  DFT  calculations  showed  that  the  dissociative 
adsorption  of  H2  on  YSZ  requires  an  activation  energy  of  at  least 
200  kj  mol-1 .  Hence  we  conclude,  in  mutual  agreement  with  the 
DFT  results  reported  in  Ref.  [23],  that  in  contrast  to  the  interaction 
of  H2  with  a  YSZ  +  O  surface  the  formation  of  OH  surface  species  as 
well  as  any  subsequent  formation  of  water  due  to  the  interaction 
of  H2  with  YSZ  is  energetically  highly  unfavorable. 

4.  Experimental  TPD  spectra  and  simulations 

In  Fig.  5a  “time  series”  of  H20  TPD  spectra  obtained  for  a 
YSZ  sample  exposed  to  ambient  air  humidity  prior  to  the  mea- 
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Fig.  3.  Illustration  of  two  low  energy  surface  reaction  pathways  I  and  II  as  obtained  in  the  present  DFT  studies  of  the  heterogeneous  H2  oxidation  reaction  on  the  YSZ  +  O 
surface.  Numbering  of  the  oxygen  surface  atoms  initially  present  at  the  YSZ  +  O  surface  is  the  same  as  in  Fig.  2  (a).  Further  details  are  given  in  the  text. 


surements  is  depicted.  The  three  consecutively  recorded  spectra 
(curves  (a)-(c))  revealed  the  decreasing  amount  of  desorbing  H20 
molecules  as  manifested  by  curve  (c),  which  was  obtained  in 
run  three,  where  practically  no  H20  desorption  signal  could  be 
observed  anymore.  However,  exposing  the  sample  (after  curve 

(c)  had  been  recorded  and  the  sample  had  been  cooled  back 
to  ca.  320 1<)  to  a  water  pressure  of  1  xl0-6mbar  for  25  min 
(corresponding  to  an  exposure  of  1100  Langmuir)  resulted  in  a 
clearly  observable  H20  spectrum  in  the  subsequent  TPD  run  (curve 

(d)  in  Fig.  5).  By  comparing  the  sum  of  the  integrated  areas  of 
the  two  TPD  spectra  ((a)  and  (b)  in  Fig.  5)  with  the  integrated 
area  of  a  CO  TPD  spectrum  obtained  from  a  platinum  (111) 
sample— exposed  to  a  partial  pressure  of  1 0-8  mbar  CO  until  a  sharp 
c(4  x  2)Pt(l  1  l)-CO  LEED  pattern  corresponding  to  a  CO  coverage 
of  0.5  monolayer  was  observed  [33]— a  value  x(H20/)  =  6  x  10-5 
could  be  estimated  for  the  mole  fraction  of  H20  solved  in  YSZ.  In 
the  latter  estimate  the  different  electron  impact  ionization  cross 
sections  of  CO  (to  CO+  [34])  and  H20  (to  H20+  [35])  were  taken 
into  account.  A  value  of  15  x  1014  atoms  cm-2  was  used  for  the 
Pt(  111)  surface  atom  density  along  with  a  YSZ  bulk  density  value 
of  6.1  gem-3.  Assuming  that  the  solution  of  H20  in  YSZ  leads  to 
the  formation  of  two  interstitial  protons  (hydrogen,  H,)  this  value 
corresponds  to  a  value  of  x(Hj)  =  12  x  10-5  for  the  mole  fraction 
of  solute  protons.  The  latter  value  is  approximately  five  times 
higher  than  the  value  of  x(H/)  =  2.2  x  10-5  obtained  in  Ref.  [36] 
for  YSZ  (10mol%  Y202)  at  comparable  H20  partial  pressure  of 
24  mbar. 

In  a  second  set  of  experiments  YSZ  samples  (after  having 
been  flashed  at  least  three  times  up  to  1200K  in  order  to  com¬ 


pletely  remove  the  YSZ  H20  contents)  were  exposed  to  a  H2/02 
atmosphere  (p(H2)=p(02)  =  1  x  10-6mbar)  before  H20  TPD  spec¬ 
tra  were  recorded.  The  results  of  these  reactive  TPD  measurements 
allow  for  comparison  with  the  results  obtained  in  a  numerical  sim¬ 
ulation  of  H20  formation  TPD  spectra  (vide  infra)  based  on  the 
H2/YSZ  oxidation  kinetics  data  used  in  Ref.  [19].  Finally  it  should  be 
noted  that  calibration  experiments  similar  to  the  one  performed  in 
the  H20  adsorption/desorption  studies  described  above  revealed 
that  YSZ  exposure  to  a  H2/02  atmosphere  also  leads  to  a  consid¬ 
erable  amount  of  YSZ  bulk  H20  formation  resulting  in  vales  of 
x(H2Oz)  =  3.2  x  10-5  and  x(H/)  =  6.4  x  10-5,  respectively,  which  are 
only  a  factor  of  two  lower  than  the  values  obtained  in  the  present 
water  adsorption  studies. 

Numerical  simulations  of  TPD  spectra  were  carried  out  using 
the  DETCHEM  software  package  [37],  which  allows  the  calculation 
of  temperature-  and  time-dependent  surface  coverages  and  surface 
species  fluxes  in  the  framework  of  a  zero-dimensional  surface  reac¬ 
tion  model  employing  thermodynamically  consistent  reversible 
elementary  reaction  kinetics. 

In  Fig.  6(a)  results  of  simulations  of  H20  TPD  spectra  obtained 
using  the  two  different  sets  of  H20/YSZ  adsorption/desorption 
and  surface  reaction  kinetics  data  previously  used  in  the  mod¬ 
eling  of  the  electrochemical  oxidation  data  obtained  for  H2/H20 
Ni-pattern/YSZ  model  anode  systems  in  Refs.  [16,20,21],  respec¬ 
tively.  The  dash-dotted  line  was  obtained  using  the  surface  reaction 
kinetics  data  as  given  in  Table  III  of  Ref.  [18].  The  dashed  line  was 
obtained  using  the  reaction  kinetics  data  given  in  Table  IV  of  Ref. 
[19].  In  each  simulation  the  thermodynamics  data  as  reported  in 
the  respective  articles  was  used.  For  comparison  a  simulated  TPD 
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Fig.  4.  Energetic  pathways  for  the  two  surface  reactions.  The  energy  pathway 
depicted  in  (a)  corresponds  to  the  reaction  sequence  Fig.  3  I.  The  energy  pathway 
reproduced  in  (b)  corresponds  to  the  reaction  sequence  Fig.  3  II. 

spectrum  obtained  using  refined  surface  kinetics  data  as  listed  in 
Table  I,  is  depicted  as  a  solid  curve.  All  simulations  were  performed 
for  a  H2  O  pressure  of  1 0-6  mbar  to  allow  for  direct  comparison  with 
the  experimental  spectrum  (open  symbols  in  Fig.  6(a)). 


Fig.  5.  Series  of  H20  TPD  spectra.  Spectrum  a:  First  spectrum  of  a  YSZ  sample 
exposed  to  ambient  air  humidity  prior  to  the  TPD  measurement.  Spectrum  b:  Sec¬ 
ond  spectrum  recorded  right  after  spectrum  a.  Spectrum  c:  Third  spectrum  recorded 
right  after  spectrum  b.  Spectrum  d:  Recorded  after  the  YSZ  sample  was  exposed  to 
a  water  pressure  of  1  x  10-6  mbar  for  25  min.  Further  details  are  given  in  the  text. 


Fig.  6.  (a)  Comparison  of  H20  TPD  spectra  obtained  in  numerical  simulations  based 
on  the  adsorption/desorption  and  surface  reaction  kinetics  data  reported  in  Ref. 
[18]  (dash-dotted  curve)  and  Ref.  [19]  (dashed  curve),  respectively,  with  the  corre¬ 
sponding  experimental  TPD  spectrum  (open  symbols)  obtained  for  a  water  pressure 
of  10-6  mbar.  The  solid  curve  was  obtained  using  the  refined  mechanism  listed  in 
Table  1,  which  was  derived  in  the  present  work,  (b)  Comparison  of  simulated  H20 
TPD  spectra  obtained  using  the  H2  oxidation  kinetics  data  (dashed  curve)  reported 
in  Ref.  [19]  and  the  refined  kinetics  data  ofTable  1  (solid  curve)  with  the  correspond¬ 
ing  experimental  H20  spectrum  (open  symbols)  obtained  for  a  H2/02  atmosphere 
(p(H2)=p(02)  =  1  x  10“6  mbar). 


Finally  in  Fig.  6(b)  an  experimental  FI20  TPD  spectrum  (open 
symbols)  measured  in  the  present  H2  oxidation  studies  in  the 
presence  of  a  H2/02  atmosphere  (p(H2)=p(02)  =  1  x  10-6  mbar)  is 
compared  with  the  corresponding  theoretical  spectrum  (dashed 
curve),  which  was  obtained  in  a  simulation  using  the  H2  oxidation 
surface  reaction  kinetics  data  of  Ref.  [19],  and  with  the  theoreti¬ 
cal  spectrum  (solid  curve),  which  was  obtained  using  the  surface 
kinetics  data  given  in  Table  I. 

The  adsorption/desorption  and  surface  reaction  kinetics  data  set 
summarized  in  Table  I  was  obtained  in  the  present  work  by  system¬ 
atically  adjusting  the  respective  pre-exponential  factor  values  (fc°) 
in  the  TPD  simulation  until  best  agreement  with  the  experimental 
TPD  spectra  was  achieved.  The  values  for  the  respective  Arrhenius 
activation  energies  (Ea)  were  taken  from  the  DFT  calculation  results 
reproduced  in  Fig.  4(a).  The  major  conceptual  difference  compared 
to  the  mechanism  of  Ref.  [19]  is  that  the  present  reaction  mech- 
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Table  I 

Summary  of  adsorption/desorption  and  surface  reaction  kinetics  parameters  derived  from  best  fit  simulations  of  the  present  experimental  TPD  curves. 


Reaction 

/<°/SM 

Ea 

H20(g)  +  []ysz  ->  H20(ad) 

(1) 

1.2  x  10~4(*,a) 

0 

H20(ad)^  H20(g)  +  []YSZ 

(2) 

8.5  x  1013  s-1 

41.7 

H20(ad)  +  [0]ysz  ->  OH(ad)  +  OH(ad) 

(3) 

1.6  x  1025  cm2  (mols)-1 

164 

OH(ad)  +  OH(ad)  H20(ad)  +  [0]Ysz 

(4) 

1.1  x  1017  cm2  (mols)-1 

173.6 

OH(ad)  +  OH(ad)  -►  H2(g)  +  [0]Ysz  +  [0]ysz(5) 

1.1  x  1010  cm2  (mols)-1 

514.4 

H2(g)  +  [0]ysz  +  [0]ysz  -►  OH(ad)  +  OH(ad)(6) 

1.2  x  10-2(*’b) 

0 

[]ysz:  YSZ  oxygen  surface  vacancy;  [0]Ysz:  “oxygen  filled”  YSZ  surface  vacancy;  £a:  activation  energy  given  in  kj  mol-1 ;  k°\  pre-exponential  factor;  So:  sticking  coefficients. 
YSZ  surface  site  density  is  1.3  x  10-9  mol  cm-2. 

A  sticking  coefficient. 
a  From  Ref.  [38]. 
b  From  Ref.  [19]. 


anism  allows  -  in  agreement  with  the  DFT  calculation  results  - 
for  the  existence  of  surface  chemisorbed  water,  H20(ad),  which  is 
formed  via  reaction  (1 )  and  which  can  subsequently  dissociate  via 
a  thermally  activated  process  (reaction  (3))  into  surface  hydroxyl 
radicals,  OH(ad). 

5.  Summary  and  conclusions 

The  comparison  between  the  H20  TPD  spectra  obtained  in 
numerical  simulations  employing  the  two  different  H20/YSZ 
adsorption/desorption  and  surface  reaction  kinetics  data  sets  of  Ref. 
[18]  (dash-dotted  curve  in  Fig.  6(a))  and  Ref.  [19]  (dashed  curve  in 
Fig.  6(a)),  respectively,  clearly  shows  that  the  kinetics  data  of  Good¬ 
win  et  al.  [19]  leads  to  a  considerable  better  agreement  with  the 
experimental  TPD  spectrum  (open  symbols  in  Fig.  6(a))  than  the 
kinetics  data  employed  by  Vogler  et  al.  [18].  In  addition,  as  can  be 
seen  in  Fig.  6(b),  the  simulation  based  on  the  FI2  oxidation  kinet¬ 
ics  data  of  Ref.  [19]  (dashed  curve  in  Fig.  6(b))  is  also  in  reasonable 
agreement  with  the  present  experimental  results  (open  symbols  in 
Fig.  6(b)).  The  good  agreement  achieved  with  the  data  set  reported 
in  Ref.  [19]  is  not  completely  unexpected,  as  it  has  been  derived 
from  experimental  measurements  of  Raz  et  al.  [15],  which  -  apart 
from  the  present  study  -  is  to  the  best  of  our  knowledge  the  only 
experimental  work  so  far,  which  specifically  has  focused  on  the 
elucidation  of  the  water  adsorption/chemisorption  and  desorption 
kinetics  on  YSZ.  Although,  the  thermogravimetry  (TG)  data  pre¬ 
sented  in  the  latter  work  is  largely  dominated  by  contributions  due 
to  physisorbed  water,  a  detailed  numerical  simulations  of  the  TG 
curves  reported  in  Fig.  17  of  Ref.  [15],  which  were  performed  by 
Goodwin  et  al.  [19],  resulted  in  a  reaction  kinetics  data  set,  which 
obviously  describes  the  actual  water  formation  and  chemisorp¬ 
tions  dynamics  on  YSZ  reasonably  well.  With  a  slight  adjustment  of 
the  kinetics  parameters  -  resulting  in  the  refined  surface  reaction 
kinetics  data  set  given  in  Table  I  -  almost  perfect  agreement  with 
the  experimental  TPD  data  could  be  obtained. 

Flowever,  as  has  been  demonstrated  in  the  present  work  as  well 
as  by  the  earlier  work  of  Wagner  [36]  the  dissociative  adsorp¬ 
tion  of  water  on  YSZ  results  in  a  noticeable  H20  dissolution  in 
the  YSZ  bulk,  either  as  interstitial  protons  (hydrogen)  or  intersti¬ 
tial  OH,  the  dynamics  of  which  has  not  yet  been  included  into  the 
detailed  chemical  mechanisms  employed  in  computational  stud¬ 
ies  of  the  electrochemical  H2  oxidation  on  Ni-pattern/YSZ  model 
anodes  although  details  about  an  estimate  of  the  current  den¬ 
sity  an  interstitial  proton  transport  mechanism  through  the  YSZ 
bulk  might  support  was  given  by  Mogensen  et  al.  [12].  In  addi¬ 
tion,  the  present  experimental  results  on  H2  oxidation  and  H20 
chemisorptions  and  dissociation  along  with  the  results  of  the  DFT 
studies  of  H2  interaction  with  YSZ  and  YSZ  +  O  model  surfaces 
clearly  demonstrate  the  important  role  the  chemical  state  of  the 
electrolyte  surface  (oxygen  poor  versus  oxygen  rich)  plays  both 
in  the  heterogeneous  H2  oxidation  and  H20  dissociation  kinetics. 


In  agreement  with  previous  DFT  studies  [23]  the  present  results 
emphasize  the  influence  the  external  oxygen  atom  supply  -  either 
by  dissociative  adsorption  of  gaseous  02  or  via  bulk  oxygen  atoms  - 
has  on  the  surface  reactivity.  Under  conditions  where  external  oxy¬ 
gen  supply  maintains  the  surface  in  an  oxygen  rich  (YSZ  +  O)  state 
the  electrolyte  surface  will  be  largely  hydroxylated  with  the  direct 
formation  of  H20  being  the  energetically  as  well  as  thermodynam¬ 
ically  favored  H2  fuel  oxidation  process.  Under  these  conditions 
the  interaction  of  gaseous  H20  with  the  electrolyte  will  result  if 
the  temperature  is  low  enough  -  as  previously  proposed  by  Raz  et 
al.  [15]  -  in  the  formation  of  molecular  water  physisorbed  on  the 
hydroxylated  electrolyte  surface.  Under  surface  oxygen  depleting 
operating  conditions,  on  the  other  hand,  the  resulting  YSZ  surface 
is  virtually  inert  toward  H2  oxidation  and  water  adsorption  fol¬ 
lowed  by  dissociation  will  become  the  dominating  process.  Hence, 
it  seems  not  unreasonable  to  anticipate  that  the  electrode  polariza¬ 
tion  (anodic  versus  cathodic)  could  have  a  pronounced  influence 
on  the  actual  surface  reactivity  of  the  electrolyte.  A  fact  that  in 
conclusion  points  towards  the  need  of  developing  more  flexible 
and  advanced  mathematical  sub-models  for  the  description  of  such 
polarization  induced  surface  reactivity  changes,  which,  after  being 
included  into  appropriate  reaction  diffusion  mechanisms  might 
contribute  finally  to  a  more  detailed  understanding  and  to  a  con¬ 
sistent  theoretical  description  of  not  only  some  selected  data  but 
all  reliable  electrochemical  experimental  data  [8]. 
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